Coherence lifetimes of excitations in an atomic condensate due to the thin spectrum 
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We study the quantum coherence properties of a finite sized atomic condensate using a toy-model 
and the thin spectrum model formalism. The decoherence time for a condensate in the ground state, 
nominally taken as a variational symmetry breaking state, is investigated for both zero and finite 
temperatures. We also consider the lifetimes for Bogoliubov quasi-particle excitations, and contrast 
them to the observability window determined by the ground state coherence time. The lifetimes are 
shown to exhibit a general characteristic dependence on the temperature, determined by the thin 
spectrum accompanying the spontaneous symmetry breaking ground state. 



PACS numbers: 03.75.Gg, 03.75.Kk, 42.50.Dv, 03.67.-a 



I. INTRODUCTION 

Although the first observation of superfluid behavior 
in helium-4 dates back to 1937 [l|, it was not until 1995 
that superfluid associated with Bose-Einstein condensa- 
tion (BEC) @ was discovered in dilute atomic gases 
While the atom-atom interactions are too complicated 
to handle in liquid state helium-4, a dilute atomic gas 
opens the possibility to construct a microscopic theory 
for superfluidity. Since 1995, atomic quantum gases have 
served as excellent fertile ground for studying quantum 
coherence properties of matter and for testing interesting 
many body theories. 

The theoretical and experimental studies of atomic 
condensates have also focused on their quantum coher- 
ence properties. Shortly after the initial discovery of 
BEC, it was understood that a finite sized condensate, in 
addition to the usual decoherence due to imperfect isola- 
tion from the environment, suffers from quantum phase 
diffusion 0, Hj] , an interaction driven decoherence due to 
atomic number fluctuations from within the condensate 
@. This study suggests a third source of decoherence, 
which we show limits the lifetime of a quasi-particle ex- 
citation from a condensate, based on the mechanism of 
thin spectrum as recently proposed and applied to any 
quantum system with a spontaneously broken symme- 
try 0, H, Q. Our work therefore constitutes a natu- 
ral application of the thin spectrum formalism to the 
highly successful mean field theory for atomic conden- 
sates, where the condensate is treated as a U(l) gauge 
symmetry breaking field. 

This paper is organized as follows: we begin with a 
review of a toy model calculation for the lifetime of the 
coherent condensate ground state as well as a squeezed 
ground state and a thermal coherent state. We then re- 
view the concept of thin spectrum and show how it is 
connected with spontaneous symmetry breaking and de- 
coherence of an atomic condensate. In sec. IV, we show 
how the quasi-particle excitations of an atomic Bose- 
Einstein condensate are affected by the thin spectrum 



associated with the ground state condensate. Finally, 
in sec. V we show how to generalize the idea of thin 
spectrum to systems with multiple broken symmetries. 
Concluding remarks are provided in sec. VI. 



II. A TOY MODEL 

The basic idea of dephasing from the ground state 
phase collapse can be understood based on the 'zero- 
mode' dynamics of a toy model [l(| EJ- The ground 
state of an N boson system is with all N bosons in the 
lowest energy eigenstate, the zero momentum state for a 
homogeneous gas. However, it cannot simply be a Fock 
state since Bose-Einstein condensation entails a definite 
phase from the broken phase U(l) symmetry, while a 
number state has no definite phase. A reasonable ap- 
proximation is to consider a coherent state occupation 
for the zero-mode with an amplitude a — \[N . Taking 
a real is equivalent to explicitly picking a phase of the 
U(l) symmetry. Because such a coherent state is not an 
energy eigenstate, it suffers phase collapse 0,0]. In this 
section, we discuss the dynamics of the associated phase 
collapse based on a simple toy model to calculate the rate 
of this collapse and study the modifications arising from 
a squeezed ground state. 



A. The lifetime for the coherent ground state 

We discuss the zero mode due to BEC, which results 
in the breaking the V(Y) gauge symmetric Hamiltonian 



hi = —a' a' aa - 
2 



(1) 



where a denotes the atomic annihilation operator for the 
condensate (zero) mode, u scales as u = uq/V with V the 
quantization volume and uq the effective interaction con- 
stant defined as uq = Aira s h 2 /M . a s is the s-wave scat- 
tering length and AI is the atomic mass, /j, is the chemi- 
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cal potential, a Lagrange multiplier for fixing the density 
of the average number of condensed particles N in the 
quantization volume V. We consider a variational, sym- 
metry breaking ground state, a coherent state satisfying 
a\z) = z\z). Such a state can be formally generated by 
the displacement operator D{z) — expfza — z*a*) acting 
on the vacuum state, or D(z)\0) — \z). Minimization of 
the mean free energy (z\Ti\z) then fixes \z\ ~ y/~N. This 
coherent state can be expanded in terms of the eigenvec- 
tors of the Hamiltonian, e.g. the Fock number states \n) 
so that 



-|z| 2 /2 



E4, 



(2) 



In this case the order parameter for BEC, is the expec- 
tation value of the annihilation operator. In the Heiscn- 
berg picture, the operator a(t) is 



a(t) = e^ m ae~ 



-Ht 



(3) 



In terms of the eigenenergy E n = ^ (n 2 — n) — /in, defined 
through 7i\n) — E n \n), for the n-th Fock state |n), one 
can easily calculate 



N exp 



(4) 



whose short time behavior is found to be 

(z\a\z) = ViVe^'e^V^' 2 , (5) 

i.e., revealing an exponential decay [1,0]. At longer time 
scale, it turns out that (z|a|z) revives due to the discrete, 
and thus periodic, nature of the exact time evolution ([4]). 

The short time decay defines a collapse-time propor- 
tional to t c ~ h/^f~Nu. The ratio of the revival time t r 
required for the order parameter scales as t r /t c = \/N, 
and becomes infinite in the thermodynamic limit. In or- 
der to get an estimate of this t c , we introduce a char- 
acteristic length scale for the harmonic trap potential as 
= y / 'h/(Mujt I ), in terms of the harmonic trap fre- 
quency wtr- Denoting the density of condensed atom 
numbers in the quantization volume as p — N/V, we 
find 



t c = 



N 



(6) 



47TiVeff W tr ' 

where we have defined N c s = pa^ a s . Assuming a typical 
situation of current experiments with N ~ 10 6 , a s = 10 
nm, <2ho = 1 M m i an d P = 10 21 m -3 , we get t c ~ 10/wtr- 
For a magnetic trap with oj tr = 100 Hz, this amounts 
to t c ~ 10 _1 seconds, clearly within the regime to be 
confirmed and studied experimentally [6| . 



B. A squeezed ground state 

The unitary squeezing operator [12] for a single bosonic 
mode is defined as 



The squeezed coherent state \a, 7) = D(a)S('f )\vac) is 
also a minimum uncertainty state, although its fluctua- 
tions in the two orthogonal quadratures are not generally 
equal to each other. Fluctuations of one quadrature are 
reduced or squeezed at the expense of the other. The 
arguments of 7 and a determine which quadrature is 
squeezed. In particular, if both 7 and a are real, then the 
state is a number squeezed state, with the uncertainty in 
atom number reduced at the cost of higher uncertainty 
in the conjugate phase variable. We expect such a state 
to have a longer life time, since the phase collapse speed 
is generally proportional to AN, which is smaller in this 
case, as have recently observed experimentally [l3|, [l4[ • 
A wide phase distribution, on the other hand, makes the 
squeezed state more similar to a Fock state which has a 
uniform phase distribution, and is less influenced by the 
decoherence effect due to the U(l) symmetry breaking 
field because of the reduced number fluctuations. 

In order to understand the essence of the above dis- 
cussion, we choose to follow similar arguments as with 
the coherent state considered previously. We will study 
the time evolution of the single mode state \a, 7) subject 
to the same U(l) gauge symmetric Hamiltonian fTJ). For 
notational convenience we define 



C = 7 



tanh(|7|) 
l7l ' 



(8) 



The Fock state expansion of the squeezed state in terms 
of this new variable is UM 



I a, 7> = ^2 A n (a, C)|j 

n=0 

= (l-ICI 2 )" 4 * 



71=0 



2C 



(9) 



(7) 



where H n is the n-th order Hermite polynomial. In the 
limit x — > co, H n (x) behaves like 2 n x n . Hence, the 
squeezed state approaches a coherent state when £ — > 0. 
The corresponding expectation value for a(t) now takes 
the form 

00 

(a l7 |a(t)ta7) =^V^TlA;A„ +1 e*^-^+ 1 )*,(10) 

where the complex nature of A n (a, makes the analytic 
evaluation of this expression nontrivial. We therefore re- 
sort to numerical studies. In a recent paper, number 
squeezing of the initial state by a factor of 10 was re- 
ported This corresponds to 7 = In 10 or £ ~ 0.98. 
In our numerical calculations, we consider the time evo- 
lution of (JTDP for a = 10 at C = 0.5 and 0.9. The results 
in Fig. [T] manifest that squeezing in the particle number 
fluctuations improves the coherence time for the conden- 
sate. The phase space distributions of the initial states 
used in Fig. [1] are displayed in Fig. [5] The longest lived 
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preparation of the condensate is the one with strongest 
squeezing in the particle number or the one with largest 
phase fluctuations. In order to examine how the phase 
distribution evolves in time, we can look the propaga- 
tion of the Q-function. We find that all coherent prepa- 
rations of the condensate eventually lose the imprinted 
phase information, and the system recovers its uniform 
phase distribution as in a Fock state. A typical result 
of our simulations is presented in Fig. [3l By initially 
preparing the condensate in a number squeezed coherent 
state, with already broad phase distribution, longer life 
times of the condensate are achieved. 




FIG. 1: The comparison of the short time decay character 
for a coherent state condensate with that for a squeezed state 
at £ = 0.5 and £ = 0.9. The parameters used are a a — 10 
nm, aho = 1 /im, n = 10 21 m~ 3 , but now for N — 100. 
In this case, the dimensionless time in units of h/u becomes 
h/u — lo^. 1 . The fastest decay (solid line) denotes the result 
for a coherent state, while the dashed (dotted) line refers to 
that of a squeezed state with ( — 0.5 (( = 0.9). As expected, 
the choice of a squeezed state with real parameters a and ( 
improves the coherence time. 



More generally, our discussions can reach beyond the 
choices of real parameters a and £. Consequently, differ- 
ent results may be expected, as we illustrate the compar- 
isons between a coherent state and squeezed states with 
C = 0.5, C = 0.5z, and C = -0.5 in Fig. g] We sec that 
the last two choices of the squeezing parameters lead to 
reduced coherence times, a result that again can be rea- 
sonably understood in terms of the increased uncertainty 
in the atom number, as it causes faster collapse. 
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FIG. 2: (Color online) The phase space distributions of the 
initial states used in Fig. [T] Curves (a), (b) and (c) corre- 
spond to ( — 0.9, £ = 0.5 and £ = 0, respectively. Although 
all of them should be centered at a — 10, they are shifted for 
convenience. 
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FIG. 3: (Color online) Time evolution of the Q-function for a 
squeezed-coherent state with a = 10 and £ = 0.5 for different 
values of tujtr- Figures (a), (b), (c) and (d) shows the Q- 
function distributions for tu) tr = 0, tw% r = 0.02, tuj tr = 0.10 
and tujtr = 0.40. It is seen that as the order parameter decays, 
the broken phase symmetry is restored, since the Q-function 
distribution becomes rotationally symmetric. 



C. A thermal coherent state 

To extend the above discussions to finite temperature 
systems, we will now introduce the thermal coherent 
state, which possesses both a thermal character as well 
as a phase. Consider the following density matrix for a 
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FIG. 4: Decay of the order parameter for the coherent state 
and squeezed states of £ = 0.5, £ = 0.5i and £ = —0.5 as 
a function of tuitr- The solid line is the coherent state, the 
dashed line is the squeezed state with £ = 0.5, and the dotted 
ones are the squeezed states with £ = 0.5i and £ = —0.5. 
Only the state with real squeezing parameter has a longer life 
time than the coherent state. 



thermal state 



ptb 



-an 



E« 



n)(n\, 



(11) 



where /3 = 1 /ksT, ks is the Boltzmann constant and T is 
the temperature. In this state pip, 7i is the Hamiltonian 
operator, H = ua 1i a^aa/2. E n is redefined, correspond- 
ing to H\n) = E n \n). p t h is a mixed state that has a 
thermal character but not a definite phase. In order to 
introduce a coherent component, and also to change the 
mean number of atoms, we can make use of the displace- 
ment operator 



= D{a)p th D\a). 



(12) 



We shall call this state a thermal coherent state, whose 
properties can be conveniently studied with the aid of 
the generalized coherent or displaced number states [l6| 



D(a)\n) — \n,a) 

oo 



m -„ Lr „ (lal 2 )ITO) 



n=0 

oo 



22 C m (n,a)\m), 



(13) 
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FIG. 5: (Color online) The phase space distributions of the 
initial states used in Fig. [4] Curves (a), (b), (c) and (d) corre- 
spond to £ = 0.5, £ = 0, £ = —0.5 and £ = 0.5i, respectively. 
Although all of them should be centered at a — 10, they are 
shifted for convenience. 



where L l k is the generalized Laguerre Polynomial. The 
thermal coherent density matrix now becomes 

oo 

p = ^e- /3B "C(a)|n)(n|^ t (a) 

n=0 
oo 

= 5>-^|n,a)(n,a| 

n=0 

= e~P E "C m (n,a)C* m ,(n,a)\m)(m'\, (14) 

nmm' 

with which we can again consider the time evolution of 
the expectation value of a(t), 

(6(t)) = ]T e- f3E "C m (n,a)Cl,(n,a)(k\m) 



nmm'k 

(m'\e^ nt ae-^ m \k), 



(15) 



n=0 



calculated according to (a(i)) = Tr(pa(t)). In the end, 
we find 

(a(t)) = ^c^ £ "C7 m+1 (n,a)C7^(n,a) 

nm 

In general, the phase factors will interfere destructively in 
the above. The thermal distribution weight e~@ En term 
determines how many different terms contribute. This 
implies that the temperature definitely leads to a reduced 
coherence time for the state. For the initial preparations 
of a condensate of 100 atoms, as depicted in Fig. El Fig. 
[S] illustrates the decay of these condensates at various 
temperatures. 
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FIG. 6: (Color online) The short time decays for thermal 
coherent states. The lines correspond respectively to T = 
1000 nK, 100 nK, 10 nK, 1 nK, and 0.001 nK from left to 
right. The humps are entirely due to the ground degeneracy 
Eo — Ej_. Even as the temperature approaches zero, the state 
(|12[) does not approach the ordinary coherent state D(a) j0). 
Instead, it is a superposition state D(a)(\0) + \l))/y/2. 



III. THE THIN SPECTRUM FORMALISM 

A. A simple theory 

By the thin spectrum, we typically refer to a group 
of states whose energy spacings are so low that they are 
not exactly controllable in any experiments. The effect of 
such states on the partition function and on the decoher- 
ence has been studied extensively, for instance see Ref. 
Q. In many-body systems, models of thin spectra arise 
quite often whenever there exists a spectrum with level 
spacing inversely proportional to the system size. These 
states with vanishing energy difference in the thermody- 
namic limit are usually beyond experimental reach and 
therefore constitute a thin spectrum. 

We begin by reviewing the ideas developed in [3] , which 
use two quantum numbers: n and to, to denote the thin 
spectrum and ordinary states. When the initial state is 
prepared at m = 0, the thin spectrum distribution will 
be a thermal one. This leads to the initial state for the 
system being 

P (t = 0) = Z- 1 ]T e-^o n) | , n ) (0, n\, (17) 

n 

where Tl\m,n) = Em \m,n) . Z is the partition func- 
tion, Z — exp (— f3 Eq ). A transformation |0,n) — > 



J2 m C m \m,n), leads to the state 

P = Z- Y e- 0E ° n) C m C m ,\m,n){m',n\, (18) 

nmm' 

which becomes 

P (t>o) = £ ^- c -^-^ 

nmm' 

C m C m ,\m,n)(m',n\, (19) 

after time evolution. When it is observed, the details 
of this density matrix cannot be seen, since the thin 
spectrum is assumed to be beyond experimental reach. 
Therefore, only the reduced density matrix, which is ob- 
tained by taking the trace of p over the thin spectrum 
states, is observed. Following the work of Q, we define 
the thin spectrum state |j'thm) by (jthin|w, n) = 6j in \m) 
where \m) denotes the ordinary observable state of a sys- 
tem. This then allows us to compute the reduced system 
state 

p (l ' Cd) = E^hm|p(t>0)bthin) 

3 

-0E M 

= E ^^e-i^-^^^CImXm'l. (20) 

mm' n 

While the diagonal elements Pram — C7 m | ^ experience no 
time evolution, the off diagonal elements suffer a phase 
collapse unless Em — E^) is independent of n. For a 
two state system (m = 0, 1), the off-diagonal element will 
decay at a rate A£ thin /£: thin with AE thin = E[ n > - E { Q n) 
and £ thin = E { n) Q. 

B. Continuous symmetry breaking and the 
Goldstone theorem 

The Nambu- Goldstone Theorem [I?) dictates the exis- 
tence of a gapless mode whenever a continuous symmetry 
is broken spontaneously. For a ferromagnetic material, 
this mode is the long wavelength spin waves For 
a crystalline structure, when the translational symmetry 
is broken, the Nambu-Goldstone mode (NGM) describes 
the overall motion of the crystal [7j . For an atomic con- 
densate, where the BEC leads to the breaking of the 
gauge symmetry, the corresponding gapless mode induces 
phase displacement of the condensate 0, [l^] . 

Consider a diagonal Hamiltonian, which may corre- 
spond to normal mode excitations with different w's, 

n = J2^kb{b kl (21) 

k 

where bk is the annihilation operator for the fc-th mode. 
As usual, the bosonic commutation relations are as- 
sumed, [bk',b\] = 5 kik > and [b k ,b k >] = [b\,b\,] = 0. If 
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there is a broken symmetry, motion along the axis of this 
symmetry will experience no restoring force, and hence 
the Hamiltonian of this mode will have the form p 2 /2I 
rather than a' a, where p is the corresponding momentum 
operator and I is the corresponding inertia mass. Hence, 
the Hamiltonian becomes 



n 



1 

2I 1 



k 



frw k b\b k . 



(22) 



The Hamiltonians for both a crystal [7j and a conden- 
sate [j| can be shown to take this form. In both cases, 
the inertia mass parameter I depends on the total atom 
number N and can either diverge or vanish in the ther- 
modynamic limit when N — > oo. 

The relationship between the Nambu-Goldstone The- 
orem and the thin spectrum is that the NGM guarantees 
the existence of a gapless mode, with the corresponding 
momentum p taking an arbitrarily small value. There- 
fore, the value ofp is always capable of giving rise to ther- 
mal fluctuations below the experimental precision and 
every NGM corresponds to a thin spectrum Q • 



C. An explicit calculation 



The Hamiltonian (|22|) is very common, Therefore, it 
is useful and instructive to calculate its time of collapse 
explicitly. According to the thin spectrum theory, the 
general state of a system takes the form \p, {Nk}) de- 
noted by two sets of quantum numbers p and {Nk}. For 
simplicity, we assume that both p and k are one dimen- 
sional quantities. Furthermore, only two different states 
of the system are considered in order to use it as a qubit. 
Assume that the elementary excitation which brings the 
system from {Nk} to {N' k } has a corresponding energy e. 
In general, such an excitation may also change the inertia 
mass / of the p term. For example, an interstitial excita- 
tion changes the total mass of the crystal 0]. Similarly, 
an excitation inside an atomic condensate can change its 
peak density, which determines the inertia mass factor in 
front of the phase coordinate [!, . Such a change is 
necessary for our mechanism of phase diffusion to occur. 
When this change to the effective mass from I to 7(1 + 6) 
is small against the small change S of the parameter, the 
off-diagonal element in equation (|2"0"1) evolves in time as 



(red) _ y- 
Pod - Z 



where E^ p) = p 2 /2I and E\ 



CiC* , (23) 



(p) 



+ p 2 /2I(l + 5). Upon 



substituting into the above, we find 



(red) _ „ 
Pod - Z 



L V 



V 21 2h I L )r 



CiC *. (24) 



Since p is continuous, its summation becomes an integral, 
or 



(red) 



Pod 



which gives 



1 



21 A hi b 



-.&(%, (25) 



|pi r d ed) | 2 = (const.)- 



Thus, the off diagonal term decays in a time 
t c — h/k B T5, 

as seen in Fig. [7] 



(26) 



(27) 




FIG. 7: The decay of 



Ip^f | as a function of t/t c . 



To apply the above result to an atomic condensate, we 
consider the relevant temperature scale at T ~ lOOnA" 
and assume that a particular observable excitation has 
5 ~ 10 _1 . In this case, we see that t c ~ 10~ 3 seconds, less 
than the life times of many observed ground states. We 
can also try to obtain an approximation to the coherence 
time of the condensate ground state. Taking the atom 
number as N ~ 10 6 , if the ground state is assumed a co- 
herent state, than the number fluctuations is of the order 
of AiV = yN. The inertia parameter I is proportional 
to I - N 2 / 5 d [H in the Thomas-Fermi limit, which 
gives S = [{N + AN) 2 / 5 - N 2 / b ] /TV 2 / 5 = 2AN/5N, or 



10 . Substituting this in, we find t c 



10" 



scc- 



onds, much larger than for the excited state, as to be 
expected. Furthermore, the result for the ground state 
life time is in agreement with our previous calculation in 
sect. II. 
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IV. QUASI-PARTICLES IN A CONDENSATE 



Substituting /xq — u oPo — uoPo/^Nq, we get 



A. A thermal state 

We now focus on the Hamiltonian of a dilute, weakly 
interacting atomic Bose gas [2l| 



H = Tn^~ n o ~ Pou n + H e 

Z1\q 



(35) 



H = J2 E ka{a k + a l+q a l- q a k a P- 



(28) 



k,p,q 



Omitting the 3 rd and 4 th order operator terms in the non- 
condensed mode (fc 7^ 0), we can partition the Hamilto- 
nian into 



For a coherent condensate state with u — UoPo, we could 
immediately get this result by neglecting the second term 
in Eq. (|34p , consistent with the non-zero momentum part 
of the Hamiltonian. 

The ground state, which we denote as |$), has Nq 
bosons in the zero momentum state and no quasi-particle 
excitations at all, i.e., for k ^ 0, 



H = H z + He, 
H z = 2^0-^0), 



(29) 



H e 



fc^O 



(E k + 2uh ) n k 



k a a 



h.c. 



or 



and 



&a|*) = 



where hk = a\a,k is the mode occupation. The above two 
parts of the Hamiltonian actually do not commute with 
each other, as we can easily check that 



|$) = 5|vac) 



(36) 



(37) 



(38) 



[H z , H e ] 



alh a 



- h. 



Neglecting the quantum nature of do in H e , we can re- 
place ri /V by po = Nq/V, and get 



Therefore, the quasi-particle vacuum state |<I>) is in fact 
a squeezed vacuum of atoms with nonzero k. 

Now, we consider a setup with n atoms in the conden- 
(30)gate mode and m quasi-particle excitations at a certain, 
single k mode, while all other modes are empty. We will 
denote such a state by \n, m) 



He 



fe#0 



efe^fe -I — (ala'_ k + h.c.) 



(31) 



no\n, m) 
n k i\n, m) 



n\n, to), 
mS k ,k'\n, to) 



(39) 
(40) 



where we have defined Ck = Ek + 2uopo — Pa with 
fio = uqPq for a coherent condensate state. In this ap- 
proximation, [H z , H e ] — at the cost of sacrificing the 
conservation of A^otal = J2k nk - 

The quadratic Hamiltonian (|31[) can be diagonalizcd 
with the Bogoliubov quasi-particles into the canonical 
form 



In the single-particle excitation regime Ek PqUq, each 
quasi-particle excitation reduces the number of conden- 
sate atoms by one. In this case, the energy of this state 
can be written as 



H\n,m) = E$\n,m) 



2(N Q - to) 



uopon + muj 



\n,m), (41) 



He 



Mo 



const., 



(32) 



with L) k = [e 2 k - ulpl] 1 / 2 [111 and b k = SakS' 1 . S is the 
multi-mode squeeze operator |23j . 

In order to conserve the particle number density in the 
condensate, we include a chemical potential term in the 
zero mode Hamiltonian 



where we simply denote u) = ujk- 

Assume the system can be initially prepared with no 
quasi-particle excitation at all, but is in a Boltzmann 
weighted distribution over the states \n, 0), i.e., 



p(t = O)ocJ2 e ~ 0E ° n) \ n ' O )( n >°\- 



(42) 



Hz = ^ <™0 ~ "0) - PoflQ- 



(33) 



The ground state of such a system will be a Fock number 
state \N ) with N = p n V/u + 1/2. We assume that 
although No and V may fluctuate, their ratio po is always 
a constant, as in the thermodynamic limit. In this case, 
H z becomes 



Hz 



uoPo 
2N Q 



{hi - h ) - p Q h Q . 



(34) 



This state will allow us to study the number fluctua- 
tions due to unknown nonzero temperature constituents 
that make up the occupations of the thin spectrum jjjj]. 
The summation index can take any positive integers and 
therefore the summation should be over < n < +00. 
However, we note that the maximum of E [ n) is at No > 
1, and because it becomes extremely small for small val- 
ues of n, we can extend the summation to be over the 
full range —00 < n < +00 and replace it with an integral 
in the continuous limit as done in the following. 
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Excitation of a quasi-particle brings each \n, 0) to 
\n, 1). The off diagonal element of the resulting state 
will evolve according to 



p od (t > 0) cx 



e (~l3u p /2No+ituopo/2hNg)n 2 +Ppou ni 

G 



exp 



,o2 2 2 

g Po"o 



2Pu po/No-2ituopo/tiNg 



v / (3u po/2N - itu p /2hN§ 



which gives 



I Pod (*) | S 



(43) 



(44) 



after omitting terms with only a phase factor. Although 
the denominator and the numerator have quite different 
forms, we find that both decay in a time proportional to 
t c ~ HNo/kBT. This is the same result that Wezel et. al. 
have found for a crystal [7j • The decay of this function 
is plotted in Fig. [8] for unit values of parameters. 



The study of temperature dependence for the damp- 
ing rates of Bogoliubov excitations of any energy has 
been carried out before using perturbation theory. A 
linear temperature dependence was found (24[, surpris- 
ingly coinciding with the linear dependence found here 
based on the decoherence of the thin spectrum. Our re- 
sult clearly would make a quantitative contribution to 
the total decay of the quasiparticles, although we note 
that our calculation is limited only to the single-particle 
excitation regime as we have used efc = Ek ^> uoPo- In 
the phonon branch corresponding to the low-lying collec- 
tive excitations out of a condensate, more complicated 
temperature dependencies may occur (25[. In contrast 
to damping mechanisms based upon excitation collision 
processes in the condensate, the thin spectrum caused 
decay rate shows no system specific dependencies, apart 
from the dependencies on temperature and the number 
of atoms. It is independent of the interatomic interaction 
strength or the scattering length, and the quasiparticle 
spectrum. This is due to the fact that thin spectrum 
emerges as a result of a global symmetry breaking in a 
quantum system so that local properties of the system 
do not contribute to the associated decay rate. 




2 4 6 8 



t/t 

c 

FIG. 8: The relative decay of the off diagonal element in equa- 
tion (|44p as a function of t/t c for unit values of parameters. 

For an atomic Bose-Einstein condensate, the relevant 
parameters are N ~ 10 6 - 10 8 and T ~ 1(T 8 - 1(T 7 K. 
These then lead to t c ~ I0 2 — 10 5 seconds, which is a 
time much larger than both the theoretical and observed 
ground state life times. However, this is the life time for a 
single quasi-particle excitation, i.e., for m = 1. It is easy 
to show that the collapse time is inversely proportional 
to to for m values not too large. An easily tractable 
excitation should have to ~ Nq and this gives t c ~ 10 -4 — 
10~ 3 seconds, much smaller than both the observed and 
expected ground state life times. 



B. A thermal coherent state 

We now generalize the above idea to a thermal coherent 
occupation of the zero-mode. The initial density matrix 
becomes in this case 

Pod (0) = Z- 1 ^-^ D{a)\n,0)(n,0\D^a) 
n 

= Z- 1 Y, e-P E o n) C m (n,a)C:An,a) 

nmm' 

|m,0)(m',0|. (45) 

The system is now brought into a superposition of no 
quasi-particle and one quasi-particle state, i.e., \n, 0) — > 
(|n,0) + |n, l))/v2- After further time evolution, the 
state becomes 

p(t) = z 2^ C m (n,a)C* n ,(n,a) 

nmm' kk f — 0,1 

e~ * (£ i m) ~ E ° m ' })t \m,k)(m',k'\, (46) 

giving rise to the reduced density matrix and its off di- 
agonal element below 

p (rcd) = z- i y E ^— \°m\ 2 

nl fcfc'=0,l 

e-H E P- E >\k)(k'\, (47) 
Pof = S e " /34 " > 1 ° l ^ ^ 1 2 e" 4 (E?) - E "" ]t - (48) 

nl 
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FIG. 9: Decay of the off diagonal element at T — 10 nK as 
a function of tuitr- Dashed line shows the decay in the case 
of the thermal coherent occupation and the solid line shows 
that in the case of the thermal occupation. 




FIG. 10: Decay of the off diagonal element for T = 100 nK 
and thermal coherent occupation of the zero mode as a func- 
tion of tuitr- Dashed line shows the decay in the case of the 
thermal coherent occupation and the solid line shows that in 
the case of the thermal occupation. 



Figs. l9l and ITOl show the early time decay at tempera- 
tures of 10 nK and 100 nK respectively. It is seen that the 
decay time for a thermal occupation, which we have stud- 
ied in the preceding subsection, exhibits stronger temper- 
ature sensitivity, whereas the decay time for the thermal 
coherent occupation is not changed very much by temper- 
ature. Therefore, we conclude that for a thermal coherent 
occupation, the main reason for the decay of the off diag- 
onal element is the decay of the zero mode distribution. 
However, if there is solely thermal occupation, no decay 



of the zero mode occurs, and the off diagonal element 
decays only because of the temperature. In the previous 
section, we have seen that the decay rate due to the thin 
spectrum and the decay rate due to the excitation col- 
lision processes show the same temperature dependence 
qualitatively. In the case of coherent thermal occupation 
of the zero mode, qualitative differences appear in the 
temperature dependence of the decay time due to the 
different decay mechanisms. Any remaining coherence in 
the zero-mode at non-zero temperatures makes the con- 
densate decay less sensitive to temperature. 



V. MORE THAN ONE BROKEN SYMMETRY 

A system may have more than one spontaneously bro- 
ken symmetry. For example, in addition to a broken 
gauge symmetry, the formation of vortices breaks the ro- 
tational symmetry of a condensate in a spherically sym- 
metric trap [26J. Furthermore, rotational symmetry can 
also be broken [27j | for a multi-component or a spinor con- 
densate (28j . When more than one continuous symmetry 
is broken, there will exist as many gapless modes as for 
the broken symmetries, each with its own thin spectrum. 
In this section, we briefly consider the effect of more than 
one thin spectrum. 

Consider a general effective Hamiltonian with two gap- 
less modes 



n 



aipf 



OLlp\ 



OL\2P\P2 



^hw k b\b k , (49) 

k 



which after a canonical transformation, reduces to 
H = aiK 2 + a' 2 p' 2 2 + ^2 tuokb\b k . 



(50) 



Without loss of generality, we use this form of the Hamil- 
tonian and henceforth omit the primes. The observ- 
able state will be denoted by n, and an easy extension 
leads to H\n,pi,p 2 ) = E ( n PuP2) \n, Pll p 2 ) with E ( ^ uP2) = 

En' * 1 + aiPi + &2P2- More generally, the primary excita- 
tion may affect both inertia terms in the two thin spectra, 
which may themselves be coupled, i.e., cti = cti(n,P2) 
and a 2 = a 2 {n,p\). Expanding around the small p\ and 
p 2 , we find around pj = 



= 4 > ) + [ ai (n,0)+ai(rc,0>2- 

[a 2 {n,0) + a' 2 {n,0)p 1 + ...]pl 
~ E^+a 1 ( n ,0)p 2 1 +a 2 (n,0)pl 



(51) 



up to the second orders in pj. Thus, we can safely ignore 
the inertia terms' dependence on other's thin excitations 
to the first approximation and let ai(n,p 2 ) — a\(n). In- 
stead of (|2"5)) we now find 



-(red) 
Pod 



= z~ 



.P1P2 



(P1)P2) j?(Pl,P2)\ 



CiC *. (52) 
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Upon substituting the approximate forms for the EjS, we 
find 

(red) x 0BCO.O) i (B (o.o) _ B (o,o> )t 

Pod —ZJ c c 

^ e -i[«i(l)-«i(0)]p?t e -|[a2(l)-«2(0)]p5t 

PI ,P2 

e -0[ai(O)pl+at2(O)pl]Q iC *^ 

p( r d ed ) = (const.)e- t ^ I) e-*/4 2) . (54) 

Thus, we see that the collapse due to different thin spec- 
tra do not influence each other severely. They combine 
to give a resulting decay with a simple single decay time 



thermal coherent ground state. The numerical results 
for a squeezed ground state reveal that phase fluctuations 
increase its coherence lifetime [IH, [li| , whereas temper- 
ature increases always decrease the lifetimes for ground 
state quantum coherence. 

The dynamics of thin spectrum are shown to lead to 
decoherence, not just on the ground state, but on quasi- 
particle excitations, or superpositions of excitations. We 
have introduced simple approximations that allowed for 
the calculations of the decoherence lifetime of the con- 
densate ground state as well as its coherence excitations. 
These calculations make possible the discussion of tem- 
perature effects in terms of the thermal and thermal co- 
herent occupations of the zero mode. We find that the 
lifetimes for these two cases are of the same order of mag- 
nitude, although the lifetime for the latter shows a weak 
sensitivity on temperature, whereas that of the former 
displays a stronger sensitivity. 



VI. CONCLUSION 

Based on a toy model calculation for the decoherence 
dynamics of a coherent ground state condensate, we have 
generalized the calculations of the dephasing times to 
cases of a squeezed coherent ground state as well as a 
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